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In this work, we have studied the thermodynamic quantities like temperature of the universe, heat 
capacity and squared speed of sound in generalized gravity theories like Brans-Dicke, Hofava-Lifshitz 
and f(R) gravities. We have considered the universe filled with dark matter and dark energy. Also 
we have considered the equation of state parameters for open, closed and flat models. We have 
observed that in all cases the equation of state behaves like quintessence. The temperature and 
heat capacity of the universe are found to decrease with the expansion of the universe in all cases. 
j^JTJ. In Brans-Dicke and f(R) gravity theories the squared speed of sound is found to exhibit increasing 

behavior for open, closed and flat models and in Hofava-Lifshitz gravity theory it is found to exhibit 
' decreasing behavior for open and closed models with the evolution of the universe. However, for flat 

universe, the squared speed of sound remains constant in Hofava-Lifshitz gravity. 
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I. INTRODUCTION 

Recently, it has become well known that the universe has not only undergone the period of early-time 
$_i | accelerated expansion (inflation) , but also is currently in the so-called late-time accelerating epoch (dark energy 
era). The unified description of inflation and dark energy is achieved by modifying the gravitational action at 
the very early Universe as well as at the very late times [1, 2]. A number of viable modified gravity theories 
has been suggested [3, 4, 5, 6, 7]. In reference [8], the connection between modified gravity and M-string theory 
' was indicated. The modified gravity gives the qualitative answers to the number of fundamental questions 
about dark energy. Indeed, the origin of dark energy may be explained by some sub-leading gravitational 
terms which become relevant with the decrease of the curvature (at late times). Moreover, there are many 
proposals to consider the gravitational terms relevant at high curvature (perhaps, due to quantum gravity 
effects) as the source of the early-time inflation. Hence, there appears the possibility to unify and to explain 
both: the inflation and late-time acceleration as the modified gravity effects [8]. Reviews on modified gravity 
are available in the references like [9] and [10]. Among the recent attempts to construct a consistent theory of 
quantum gravity, much attention has been paid to the quite remarkable Hofava-Lifshitz quantum gravity [11]. 
An extensively studied generalization of general relativity involves modifying the Einstein-Hilbert Lagrangian 
in the simplest possible way, replacing R — 2A by a more general function f(R) [10, 11, 12]. Recently the 
modified Hofava-Lifshitz f(R) gravity has been proposed in ref.[l]. Discussions on Hofava-Lifshitz gravity have 
been made in references [13, 14]. The basic idea of Hofava-Lifshitz gravity is to modify the UV behavior of 
the general theory so that the theory is perturbatively renormalizable [14]. However this modification is only 
possible on condition when we abandon Lorentz symmetry in the high energy regime [14]. In reference [15], 
the very interesting physical implications of Hofava-Lifshitz gravity are summarized as: (i) the novel solution 
subclasses, (ii) the gravitational wave production, the perturbation spectrum, (iii) the matter bounce, (iv) the 
dark energy phenomenology, (iv) the astrophysical phenomenology, and (v) the observational constraints on 
the theory. Recently scalar-tensor theories have received renewed interest. The Brans-Dicke theory [16] is the 
simplest example of a scalar-tensor theory of gravity. In Brans-Dicke theory, Newton's constant becomes a 
function of space and time, and a new parameter uj is introduced. General relativity is recovered in the limit 
ui — > oo [17]. Interacting dark energy [18, 19] and holographic dark energy [20, 21, 22, 23, 24] models have 
been considered in Brans-Dicke theory. Brans-Dicke scalar field as chameleon field has been considered in the 
references [25] and [26]. 

A profound connection between gravity and thermodynamics was first established by Jacobson [27], who first 
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showed that the Einstein gravity can be derived from the first law of thermodynamics in the Rindler spacetime. 
Thermodynamic aspects of the cosmological horizons have been reviewed in [28] and [29]. Investigating the 
generalized second law (GSL) of thermodynamics in gravity has gained immense interest in recent years. A 
plethora of papers have studied the thermodynamics in Einstein gravity theory [30, 31, 32, 33, 34, 35]. As the 
modified theory of gravity was argued to be a possible candidate to explain the accelerated expansion of our 
universe by various authors [36, 37, 38], thus it is interesting to examine the GSL in the extended gravity theories 
[39, 40, 41, 42]. Thermodynamics has been studied in the brane world scenario [43, 44, 45, 46], Hofava-Lifshitz 
gravity [47, 48, 49], Brans-Dicke gravity [50, 51, 52] and in f(R) gravity [53, 54, 55]. Extending the study of [15], 
two of the authors of the present paper, examined the validity of the GSL in various cosmological horizons of 
a universe governed by the Hofava-Lifshitz gravity and the GSL was proved to be valid in different horizons [56] . 

In the present work, we have studied the thermodynamic quantities of the universe in generalized gravity 
theories like Brans-Dicke, Hofava-Lifshitz and f(R) gravities. Instead of investigating the validity of the laws 
of thermodynamics, we have tried to investigate how the thermodynamic quantities like heat capacity (C„), 
temperature T and squared speed of sound behave during the evolution of the universe governed by the 
said gravity theories. In addition to this, the equation of state parameters have also been studied for all of the 
said gravity theories. Organization of the rest of the paper is as follows: In section II, we have discussed the 
thermodynamic quantities. In sections IIIA, IIIB and IIIC we have discussed the thermodynamic quantities 
under Brans-Dicke, Hofava-Lifshitz and f(R) gravity theories respectively. Finally, in section IV, we have 
discussed the results. 



II. GENERAL DESCRIPTION OF THERMODYNAMIC QUANTITIES 

The Einstein field equations for homogeneous, isotropic FRW universe are given by [57] (choosing c = 1) 



, k 8ttG 



(1) 



and 



H-t = -AitG{p + p) (2) 

where H(— ^) is the Hubble parameter and k = 0, — 1,+1 denote the curvature index for flat, open and 
closed universe respectively. Here, p and p denote the energy density and pressure of the universe. The energy 
momentum tensor T/f is conserved by virtue of the Bianchi identities, leading to the continuity equation [57] 

p + 3H(p + p) = (3) 
where p is the isotropic pressure and p is the energy density of the fluid defined by 



Here, U is the internal energy and V is the volume of the universe. 

We consider the FRW universe treated as a thermodynamical system. Then from Gibb's equation of thermo- 
dynamics, we have [33] 



TdS = d(pV) +pdV = d((p + p)V) - Vdp (5) 

where S is the entropy, T is the temperature and V is the volume of the universe. The integrability condition 
of thermodynamic system is given by [58] 



which leads to the relation between pressure, energy density and temperature as 



dp= P -^dT (7) 



dS = d[<P±2>?) (8) 



From (5) and (7), we get 



T 

and integrating, we can obtain the expression of the entropy as (except for an additive constant) 

However, for adiabatic process entropy is constant and consequently, the equation (5) becomes 

d[(p + p)}=Vdp (10) 

Relation (9) can also be obtained using (7) into (10). Hence for adiabatic process equation (9) may be 
considered as the temperature defining equation as 

t= (e±pK 

The square speed of sound and heat capacity are defined by 

< - 1 <-> 

and 

d 

C v — Vq^; (when entropy S is constant = So, say) (13) 

These thermodynamic quantities would be investigated for their evolution with the expansion of the universe 
in the subsequent sections. 



III. THERMODYNAMIC QUANTITIES IN GENERALIZED GRAVITY THEORIES 

A. Brans-Dicke Theory 

The Jordan-Ficrz-Brans-Dicke theory (heretofore, we will call it Brans-Dickc (BD) theory for simplicity) is 
the simplest example of a scalar-tensor theory of gravity. A brief introduction of the BD theory has been 
presented in the previous section. The Lagrangian density for the Brans-Dickc theory is 



-^R+^d^ + Cr, 



(14) 
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where <f> is the Brans-Dicke field, and L m is the Lagrangian density for the matter fields. The self-interacting 
BD theory is described by the Jordan-Brans-Dicke (JBD) action (choosing c = 1) as: 



16tt 



- V(<f>) + l&ir£ r . 



(15) 



where V{4>) is the self-interacting potential for the BD scalar field 4> and u){<j>) is modified version of the BD 
coupling parameter which is a function of <f>. In this theory ^ plays the role of the gravitational constant G. 
This action also matches with the low energy string theory action for u = — 1. The matter content of the 
Universe is composed of matter fluid, so the energy-momentum tensor is given by 



T JZ = (P + P)u^u v + p 

where is the four velocity vector of the matter fluid satisfying u fl u l/ 
density and isotropic pressure. 

From the Lagrangian density we obtain the field equations 



(16) 

-1 and p, p are respectively energy 



and 



8tt 



1 u riM V W 

[<P, K u - g^v u <P\ ^p~9^ 



(17) 



U6 = 



8ttT 



1 



2V(0) 



dv{4>) 



3 + 2cj(0) 3 + 2w(0) 
where T = T™ v g^ v . Equation (17) can also be written as 



dui(tp) 
d(j) 



3 + 2u(4>) 



(18) 



(19) 



where T„„ can be treated as effective energy momentum tensor. The line element for Friedman-Robertson- 
Walker space-time is given by 



ds 2 



-dt 2 + a 2 (t) 



dr 2 



1 — kr 2 



+ r 2 (d6 2 + sm 2 6d(j) 2 ) 



(20) 



where, a(t) is the scale factor and k (= 0,-1, +1) is the curvature index describe the flat, open and closed 
model of the universe. 

We are considering the universe filled with dark energy (with energy density p£>) and dark matter (with 
energy density p m ). As we are not considering interacting situation, the conservation equations are separately 
satisfied for dark matter and dark energy. Thus 



and 



p D + 'ZH{p D +p D ) = 



p m + 3H(1 + w m )p m = 



(21) 



(22) 



Solving the conservation equation for dark matter, we get the density and pressure of dark matter as p m — 
p mo (l + z )3(i+» m ) an d pm = Pmo (i + Wm )(l + 2)3(i+t«m). Defining pi = ?f and Pl = *f the Einstein's field 
equations can be written as 



5 



-0.990 F 



-0.992 - 



-0.994 



-0.996 



-0.998 




Fig.l 

Fig. 1 shows the EOS parameter w to tai = for fc = — 1 (the red line), k = 1 (the green line) and k = (the blue 

line) for Brans-Dicke model where dark energy and dark matter satisfy the conservation equation separately. We have 

taken a = 3, B = 0.1, p m o = 0.23, w m = 0.003. 



H 2 + — = —{pi+p D ) 

Qi o 



(23) 



H 



-4tt(pi +Pi+ Pd +Pd) 



(24) 



where 



u ^ 3 R 4> [ y(0) 

167T </> 2 87T (A 167T0 



(25) 



and 



^ 2 | g 10 V^) 

167T 2 An (f) 8tT (f) 167T0 

To find the thermal quantities we use the choices of <fi and V as 



(26) 



4> = cj) a a , 



V = V t 



(27) 



Using EOS for dark energy pu = wpn and the solution of (21) pn = poo - 3 ( 1+ "'*' in the field equations we 



get 



H 2 = kAa- 2 + Bfl-^T 1 + Co- a - 3 < 1+1 °) 



(28) 



(29) 
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Fig.4 

Fig. 2 shows the heat capacity C v for k = — 1 (the thick line), k = 1 (the dotted line) and k = (the broken line) for 
Brans-Dicke model where dark energy and dark matter satisfy the conservation equation separately. We have taken 

a = 3, B = 0.1, p m o = 0.23, w m = 0.003. 
Fig. 3 shows the temperature T with evolution of the universe for k = — 1 (the red line), k = 1 (the green line) and 

k — (the blue line) for Brans-Dicke model. We have taken a = 3, B — 0.1, p m o = 0.23, ui m = 0.003. 
Fig. 4 shows the squared speed of sound vl with evolution of the universe for k = — 1 (the red line), fc = 1 (the green 
line) and k — (the blue line) for Brans-Dicke model. We have taken a = 3, B = 0.1, p m o = 0.23, w m = 0.003. 



where, 



C = — 



167r(l+m)ppo 



A 2 



A = i 

a((l+w)a-l)-(2+a) 
2+a(l-a(l+^))A . R 

2+ck ' 1 2+a 

6(f> +AB 2 /B; B 2 =6 PD0 {l + a-uja 2 /6)B; C 2 



</,o[(2+a)(-a-3(l+t«))+2a((l+cj)Q-l)l 
(l-a(l+i^))B , q _ a{l-a{l+io))<t> C-&n(l + 



>)po 



(2+a)0 o 

B 2 C/B ~ 16irp D0 



(30) 



Equation of state parameter Wtotai = is computed for fiat, open as well as closed universes and are 

plotted against redshift z in figure 1 and it is found that in all of the above cases the EOS parameters are 
staying above —1, which indicates quintessence- like behavior. The behaviour of the EOS parameters further 
indicate that the energy density is increasing with evolution of the universe irrespective of it curvature. 

We replace p and p by (p\ + pd) and (p\ + pn) in equation (11) we get the temperature T, which is used in 
equation (13) to get the heat capacity C v . Heat capacity C v is computed for open, closed and fiat universes 
and are plotted in figure 2. This figure shows that for all of the three universes, the heat capacity is decreasing 
with increase in the redshift. This means that the heat capacity is increasing with evolution of the universe 
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irrespective of its curvature. Also, we present the temperature T against redshift z for all of the curvatures. 
We find that the temperature is decreasing with evolution of the universe. 



B. Hofava-Lifshitz Gravity 

Thermodynamics in cosmology has been extensively studied either in Einsteins theory of gravity or in modified 
theories of gravity. In this section, we shall generalize such studies to the Hofava-Lifshitz (HL) Cosmology. An 
exhaustive review of HL cosmology is available in [59] . We briefly review the scenario where the cosmological 
evolution is governed by HL gravity. The dynamical variables are the lapse and shift functions, N and Ni 
respectively, and the spatial metric . In terms of these fields the full metric is written as [60] 

ds 2 = -N 2 dt 2 + g tJ (dx l + N l dt){dx 3 + N ] dt) (31) 

where indices are raised and lowered using <?y . The scaling transformation of the coordinates reads: t — » IH 
and x l — > lx l . 

The action of the HL gravity is given by [60] 



C -- 
Ci = V9 N 



I = dt J dtd 3 x(C + Ci 
^gN\^{K rj K^ -\K 2 )^ 



Cm) 

« 2 M 2 (Afl-3A 2 ) 
8(1-3A) 



(32) 



«> 2 (1-4A) 
32(1-3A) 



where, k 2 , A, /i, lu and A are constant parameters, and Cij is Cotton tensor (conserved and traceless, vanishing 
for conformally flat metrics). The first two terms in Co are the kinetic terms, others in (Cq+C\) give the potential 
of the theory in the so-called "detailed-balance" form, and C m stands for the Lagrangian of other matter field. 
Comparing the action to that of the general relativity, one can see that the speed of light and the cosmological 
Newtons constant are 
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Fig. 6 shows the behavior of temperature T with evolution of the universe for fc = — 1 (the red line), k = 1 (the green 
line) and k — (the blue line) for Hofava-Lifshitz gravity. We have taken A = 1.2, p — 1.03, w m = 0.03 and p m o = 0.23. 
Fig. 7 shows the behavior of heat capacity C v with evolution of the universe for k = — 1 (the red line), k — 1 (the green 
line) and k — (the blue line) for Hofava-Lifshitz gravity. We have taken A = 1.2, p — 1.03, w m = 0.03 and p m o = 0.23. 




Fig. 8 shows the squared speed of sound with evolution of the universe for k — — 1 (the red line), k = 1 (the green 
line) and k — (the blue line) for Hofava-Lifshitz gravity. We have taken A = 1.2, /i = 1.03, w m = 0.03 and p m o = 0.23. 



C= ~VT^3A' Gc = 16^(3A - 1) (33) 

It may be noted that when A = 1, Cq reduces to the usual Lagrangian of Einsteins general relativity. Thus, 
when A = 1, the general relativity is approximately recovered at large distances. 

As we are considering dark energy with dark matter the conservation equations are given by (21) and (22). 
The field equations are 



8ttG c 



{Pm + PD) 



(34) 
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and 



H+^H 2 + ^ = -^G c { Pm + PD ) (35) 



where 



3k V^ 2 3k V A 2 1 /3fc 2 \ 

PD = 8(3A-1K + 8(3A = 1) = T^G C [a^ +3A ) (36) 

and 

PD - 8(3A - l)a 4 8(3A - 1) ~ 16nG c U« 4 / 

Using the solution for the conservation equation for dark matter given in (22) we get the total energy density 
as a function of rcdshift (z = - — 1) as 



(38) 



P(z) = Pm + PD = p mo (l + zT {1+w ^ + — ^- ( ^^^ + 3A) 
Similarly, the total pressure as a function of redshift z is 

p(z) = Pm +Pd= Pmo w m {l + z) 3 d+^) + _L- (^ fc2(1 + z)4 - 3a) (39) 
The squared speed of sound is given as a function of z by 

2, \ = fc 2 (l + z)K 2 H 2 + 6(1 + Zf w ™{-1 + ZX)W m pm0{l + W m ) 

™ 3{fc2(l + z) K V + 2(l + z) 3 --(-l + 3A)(l + W ;„ l )p m o} 1 j 

and the heat capacity becomes 

fi r \ _ 36'o{fc 2 (l + z)«V + 2(1 + zf w ™{-l + 3A)(1 + w m )p m0 } , , 

v[Z) F(l + z)«V + 6(l + z) 3 ^(-l + 3A) Wm p m0 (l + Wm ) { ' 

The thermodynamic quantities expressed above are now plotted against redshift to see their behavior with 
the evolution of the universe. In figure 5, where we have plotted the equation of state parameter for the Hofava- 
Lifshitz gravity, we see that the behavior is like Brans-Dicke theory. It is staying above — 1. However, at lower 
redshifts the equation of state parameter is tending to —1. However, it never crosses —1. Like Brans-Dicke, 
this behavior remains the same for fiat, open and closed universes. In figure 6 we find that in the case of flat, 
open and closed universes, the temperature T is decreasing with the evolution of the universe. From figure 6 
we see that the heat capacity C v is increasing as we are approaching towards the lower redshifts. From figure 
7 we understand that for open and closed universes, the squared speed of sound v 2 decreases with the evolu- 
tion of the universe. However, for flat universe, the v 2 remains constant throughout the evolution of the universe. 



C. f(R) Gravity 

Motivated by astrophysical data which indicate that the expansion of the universe is accelerating, the modified 
theory of gravity (or f(R) gravity) which can explain the present acceleration without introducing dark energy, 
has received intense attention. Extensive review of f(R) gravity is available in [61]. The action of f(R) gravity 
is given by [62] 



^matter 



(42) 
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where g is the determinant of the metric tensor g^, C mat ter is the matter Lagrangian and k 2 = 8nG. The/(i?) 
is a non-linear function of the Ricci curvature R that incorporates corrections to the Einstcin-Hilbert action 
which is instead described by a linear function f(R). The gravitational field equations in this theory are 



H 2 + ^ = 



3f'(R) 



(p + Pc) 



(43) 



H 



2f'(R) 



(P+P + Pc+Pc 



(44) 



where p c and p c can be regarded as the energy density and pressure generated due to the difference of f(R) 
gravity from general relativity given by [61] (choosing G — 1) 



8tt/' 



f Rf ' -3Hf"R 



(45) 



1 



Pc = 



8tt/' 



- n Rf +f"R + .f"'R 2 +6f"R 



(46) 



where, the scalar tensor R = -6 (^H + 2H 2 + ^ . 

As we are considering both dark matter and dark energy, the Friedman equations take the form 



2 k _ 8-7T 
ti H 2 — —p total 

CI O 



(47) 



k 

H J = -^{ptotal +Ptotal) 



(48) 



where, 



Ptotal = Pi + Pc 



Ptotal =Pl+Pc 



(49) 



with pi — £jr and p 1 — ^p-. As there is no interaction, like the previous two cases the dark energy and dark 
matter satisfy the conservation equation separately. Therefore, we have the density and pressure of dark matter 
as p m = Pm (l + z) 3 ( 1+Wm ) and p m = p mo (l + w m )(l + Z ) 3 ( 1 + M ' m ). I n the present section, while considering the 
f(R) gravity, we have illustrated with a solution, 



A 

f(R) = /3R + aR m , R=— , m > 1, a>0, /3>0 



(50) 



Also we have chosen n(m — 1) = 1. 

Using the above form of f(R) in (44) and (45) we have computed the temperature T, squared speed of sound 
v 2 and heat capacity C v as functions of the redshift z as follows: 



Pi = 



(l + z 



,3(1+U) m ) 



PmO 



m(A(l + z) n )- 1+m a + , 



(51) 



Pl m(A(l + z) n )- l + m a + [3 



(52) 



11 




-0.5 0.0 0.5 1.0 1.5 2.0 



Fig. 11 

Fig. 9 shows the plot of heat capacity C v against redshift z in f(R) gravity. We see that C v is increasing with the 
evolution of the universe. Here a = 2.10132, /3 = 2.10101, w m — 0.003, m = 2.3, p m o = 0.23 and the red, green and 

blue lines correspond to k = — 1, 1, respectively. 
Fig. 10 shows the plot of temperature T against redshift z. We find that the temperature is decreasing with the 
evolution of the universe in f(R) gravity a = 12.1, /3 — 10.1, w m — 0.003, m = 12.3, p m o = 0.23 and the red, green 

and blue lines correspond to k = — 1, 1, respectively. 
Fig. 11 shows the plot of the squared speed of sound v'j, against redshift z in f(R) gravity. We find that is increasing 
with the evolution of the universe. Here, a = 12.1, /3 = 10.1, w m — 0.003, m = 1.9, p m o = 0.23 and the red, green and 

blue lines correspond to k = — 1, 1, respectively. 



q(-l+m)^(l + Z r f mn(2m - 4)(1 + z) 4 ' 2 " - ggg + ^T"" ) 

Pc 167r(mA- 1 + m (l + 2r)™(- 1 + m )a + /3) | + A 2 M 5 ) 



P C = {"3^ + 6(-2 + m) mn 2 (l + n) 2 (l + z) 4 -»+ 

+6mn(l + z) 2 -"(-5 +n + z + nz)^C x - ^7+ ^gjgf^ } 



(54) 



Using (51), (52), (53) and (54) we get temperature T, squared speed of sound v 2 and heat capacity C v as 
functions of redshift z in the following forms 

rj, Pl+Pl+Pc+Pc 

= (l + z)35 (55) 
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where 



= [l6(m(A(l + z) n )- 1+m a + /3) 2 x 

An(l+z)" 1 + 4 "(-3A+6(-2+m)mri 2 (l+n) 2 (l+z) 4 ""+6mri(l+z) 2 ""(-5+n+z+riz)y^Ci _ 



fc , A(l + s)-4+7i 
(l + z)^ 3(-4+n) 



24m7r 48m7r 



3rnn(l+z) 2 +"(-5+rt+z+nz)f 2fc , +k(l+^)~ 5+ ") 

-6(-2 + m)m(-4 + n)n 2 (l + n) 2 (l + z) 3+n + W)3 I >- + 

In k i ^(i+»)~ 4+ " 

Y Gl_ <m) 2 " + 3 ( -4+„) 



and 



6mn(l + n)(l + z) 2+ " - + - 6m(-2 + n)n(l + z) 1+ "(-5 + n + z + nz)x 

'(57) 



r _ a_ A(l+z)- 4 +" \ 3u' m (l+«>, 71 )(l+z) 2 + 2 '*+ 3m '"p r71 o A(-l+m)wm tt . r71 (l+z) 2 +"+ 3m "'(A(l+z)")"'ap m0 
U1 (1 + z) 2 ~T 3(-4+n) J 1 m( J 4(l+z)")- 1 +™Q+/3 (m(A(l+z)")™a+A(l+z)"/3) 2 



f3Ci 

-(-1 + m)mn(i(l + z)") m ( 1 + ^p— 5 — - 




{(-1 + m){A{\ + z) n )- 1+m a + m(A(l + z) n )- 1+m a + /3}] + 

2(-2+m)(-l+ m )mn(l+ Z ) 3 - 2 "(A(l+ Z )")™^-6(-2+n)Ci + ^ 4 ^i^ 

16(-1 + m)mn(l + z) 2+3 ^™(A(l + z) n )- 1+m a Pm o + 48(1 + w m )(l + z) 2+3tu -"(™(A(l + z)")- 1+m a + /3)p m0 

(58) 

and 

c - - 18 < 59 > 

where 

(-l+m)r r m(l+z)- 1 + "(A(l+z)")- 1 +"|A 2 +2(-2+m)mn(l+z) 4 - 2 "^3C 1 - TT j^+ A(1 + ^ tl 4+ " ^| 

= 16A7r(m(A(l+z)")- 1 + m a+ / g) aX 

f 1 «(-l+m)(.4(l+z)")- 1 + '" l 
\ ±T (m(A(l+z)")- 1 + ™a+ / 3) J "r" 

(-l +m )(A(l+zn m {2(-2+m)mn(l + z) 4 - 2 "( IT ^^ 

^ 16Tr(m(A(l+z)")- 1 + m a+/3)A 2 

A(-l+m)m«(l+z)' 1+ " +3(1+m "' ) (A(l+z)")- 2 +"'ap m o 3(l+u. r71 )(l + z)- 1 + 3 ' 1 +"'"'p m o 
(m(A(l+z)")- 1 +™a+/3) 2 ^ m(A(l+z)™)- 1 +™a+l3 

(60) 
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Fig. 12 shows the evolution of the equation of state parameter w to tai with the evolution of the universe in f(R) gravity. 
We find that for k = — 1, 1, the equation of state parameter w to tai > — 1. The indicates quintessence era. 
Here, a = 10.32, j3 = 10.01, w m = 0.003, m = 12.3, p m o = 0.23 and the red, green and blue lines correspond to 

k = — 1, 1, respectively. 



and 



Sotl+z) 7 



,-r( 



A(l+z) 2 ™ -3A+6(-2+m)mn 2 (l+n) 2 (l+z) 4 ~"+6rrm(l+z) 2 ~"(-5+n+z+riz) A / Ci 



(l+z)* 3(-4+n) 



48m7r 



I -mil \: j A 2 +2(-2+m)mn(l+z) 4 " 2 " ^3Ci-. 

167rA a (m(A(l + z)")- 1 + m a+^) 



3fc , Atl+z) 



J_ (1 + Z) / ' . . : ■ _ > : . 1 1 " - 1 

" r m( J 4(l + z)' l )- 1 + m Q+^ "r" m(A(l+z)")- 1 + '"a+^ 



(l+z) 3(1+m '" ) ' J "'P 



+ 



S (l+z) 6 



< n , >2n / „ 3mrt(l+z) 2 ~"(-5+n+z+7tz) ( 2k , + j r A(l + z'r 5+ " S ) 

A(1 + z) -6(-2 + m)m(-4 + n)n 2 (l + n) 2 (l + z) 3 "" + -, W)3 3 >- 



Ci- 



fc A(l + *)— i+" 

(I + 2) 2 ^ 3(-4+n) 



6mn(l + zf-nJCx - + ^+^^ ((1 + „)(l + z) - (-2 + n)(-5 + n + z + tw))) + 

An(l+zr 1 + 2 "^-3A+6(-2+m)mn 2 (l+n) 2 (l+z) 4 -''+6m7t(l+z 



(-l+m) mn (A(l+z) n ) m a\ A 2 +2(-2+m)mn(l+z) 4 - 2 "( 3Ci-^j^ 7 + ' 4(1 +" 4 ) + ~ 4+ " J I , 

16A 2 (l+z)7r(m(A(l+z)")- 1 +'"a+ / 3) \~ 



(-2+m)(-l+m)mri(l+z) 3 " ! "(A(l+z)°) m -6(-2+n)Ci + 



6fc(4-5n + ^)(l + z) z -.4n(l + z) r 
(-4+n)(l + z) 4 



(-l+m)(A(l+z)")- 1 + '" 
(m(A(l+z)")-i + "'Q+/3) 



, 3(l+m m )(l + z) fa 
m(A(l+z)")- 1 + '"a+^ + 



8^7r(m(A(l+z)™)- 1 + ™a+,S) 
3m m (l+«- m )(l+z) 2 + 3 """p m o A(-l+m)mn(l+ 2 ) ;i - | '" w "-"(A(l+z)")' n ap m i) /, , v 
m(A(l+z)")- 1 + m Q+/3 (m(A(l+z)") m Q+A(l+z)",3) 2 V 1 T %i 

(61) 

The thermodynamic quantities expressed above are now plotted against redshift z to sec their behavior 
with the evolution of the universe. We proper choice of the parameters we plot all of the quantities in 
figures 8, 9, and 10 respectively. In figure 8 we find the increasing behavior of the heat capacity with the 
evolution of the universe. This behavior remains the same irrespective of the curvature of the universe. 
From figure 9 we see that as the universe in evolving, the temperature T is decreasing. Here also we 
get the same behavior for open, closed and flat universes. It may be interpreted that the temperature of 
the universe decreases as it expands under f(R) gravity. In figure 8 we plot the squared speed of sound 
i> 2 . Like the temperature, v 2 s is decreasing with the expansion of the universe under f(R) gravity. The 
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choices of the parameters are mentioned in the figure captions. Behavior of the equation of state parameter 
Wtotai is observed in figure 11. Throughout the evolution of the universe Wtotai > — 1- This indicates 
quintessence like behavior of the equation of state parameter. Therefore, we see that in f(R) gravity, 
where we are considering the coexistence of dark energy and dark matter with very small pressure without 
interaction, the equation of state parameter behaves like quintessence era. This holds true for flat, closed 
as well as open universes. However, it also discerned that Wtotai is gradually increasing in the negative direction. 



IV. DISCUSSIONS 

In the present work, we have considered modified gravities as Brans-Dicke, Hofava-Lifshitz and f(R) 
gravities. Various thermodynamic quantities like temperature, heat capacity and squared speed of sound have 
been investigated for all the gravity theories. In each case we have considered that the universe is filled with 
dark matter and dark energy which arc not interacting. Prior to evaluating the thermodynamic quantities 
we have studied the behaviors of the equation of state parameters. In figure 1 we see that in the case of 
Brans-Dicke gravity theory, the equation of state parameter is staying above —1 throughout the evolution 
of the universe. This indicates that the equation of state parameter is behaving like quintessence in this 
case irrespective of the curvature of the universe. In figure 5, where we have plotted the equation of state 
parameter for the Hofava-Lifshitz gravity, we see that the behavior is like Brans-Dicke. It is staying above 

— 1. However, at lower redshifts the equation of state parameter is tending to —1. However, it never crosses 

— 1. Like Brans-Dicke, this behavior remains the same for flat, open and closed universes. Similar behavior 
of the equation of state parameter is discernible in f(R) gravity also. The evolution of the equation of state 
parameter for f(R) gravity has been presented in figure 12. From the figures 3, 6 and 10 we see that the 
temperature T of the universe is decreasing with evolution of the universe in Brans-Dicke, Hofava-Lifshitz and 
f(R) gravities respectively. This behavior remains the same in open, closed and flat universes. From figures 
2, 7 and 9 we find that the heat capacity C v of the universe increases with the evolution of the universe. 
Moreover, in all of the cases C v remains at the positive level throughout the evolution of the universe. We have 
also investigated the squared speed of sound v 2 in all of the cases. From figure 4 we see that for closed universe 
v 2 starts decreasing from redshift —0.2. However, up to z = —0.2 it has gradually increased. In the cases of 
open and flat universes, v 2 has an increasing behavior throughout the evolution of the universe. From figure 8 
we see that in Hofava-Lifshitz gravity v 2 has a decaying behavior throughout the evolution of the universe in 
the case of open and closed universes. However, for flat (k = 0) universe, the squared speed of sound remains 
constant throughout the evolution of the universe. From figure 11 we understand that v 2 increases throughout 
the evolution of the universe irrespective of the curvature of the universe. 
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